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Abstract 

In a two-dimensional toy model, motivated from five-dimensional heterotic M-theory, we study 
the collision of scalar field kinks with boundaries. By numerical simulation of the full two-dimensional 
theory, we find that the kink is always inelastically reflected with a model-independent fraction of its 
kinetic energy converted into radiation. We show that the reflection can be analytically understood 
as a fluctuation around the scalar field vacuum. This picture suggests the possibility of spontaneous 
emission of kinks from the boundary due to small perturbations in the bulk. We verify this picture 
numerically by showing that the radiation emitted from the collision of an initial single kink eventually 
leads to a bulk populated by many kinks. Consequently, processes changing the boundary charges 
are practically unavoidable in this system. We speculate that the system has a universal final state 
consisting of a stack of kinks, their number being determined by the initial energy. 
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1 Introduction 



In Ref. we presented a five-dimensional brane-world model, closely related to five-dimensional het- 
erotic M-theory [21 El HI El El , where M-theory five-branes are modelled as kink-solutions of a bulk scalar 
field theory. We were particularly interested in studying collisions of such kinks with the boundaries, a 
process which, in the original M-theory model, may lead to a topology-changing so-called small-instanton 
transition [71 [HI due to absorption of the five-brane. The analysis of Ref. 1. was mainly performed in 
a moduli space approximation. We concluded that the colliding kink was absorbed by the boundary, 
however, its final fate could not be determined due to a break-down of the moduli space approach. 

In this paper, we will focus on an even simpler, two-dimensional model which captures the essential 
features of its five-dimensional cousin. The model consists of a (l + l)-dimensional scalar field theory for 
a single real scalar and an associated potential V with (at least) two minima to allow for kink solutions. 
The spatial direction is taken to be a line segment with the boundary conditions being provided by the 
"superpotential" W associated to V. To simplify the problem we have not included gravity and 
any other gravity- like fields which were present in the five-dimensional model pp. 

Our main goal is to determine the final outcome of the kink-boundary collision using this simple two- 
dimensional model. In particular, we would like to clarify whether the kink is indeed absorbed or, rather, 
reflected by the boundary. While the former process would lead to a change in the boundary charge (as 
measured by the superpotential value) and the number of kinks present in the bulk, the latter process 
would conserve those numbers. We will mainly approach this problem as one posed within relatively 
simple two-dimensional scalar field theories whose basic properties we are investigating, although the 
possible relation to M-theory and topology-changing phenomena is in the back of our minds. 

The plan of the paper is as follows. In the next Section, we will set up the model and derive its 
basic kink solution. Section 3 reviews the moduli space description of kink evolution and kink-boundary 
collision. The perturbation spectrum on the kink background is analysed in Section 4. In Section 5, we 
present a full numerical analysis of the two-dimensional model to determine the final outcome of the 
collision process. An analytic interpretation of the numerical results is given in Section 6. In Section 
7, we present a modified moduli space picture which incorporates some, although not all, features of 
the collision process. We conclude in Section 8 by summarising and presenting some results about the 
long-time evolution of the system. 

2 The model 

In this section, we will set up the 1 + 1-dimensional model which will be the basis of this paper. This 
model is designed to capture the essential features of five-dimensional heterotic M-theory [2] and the 
related defect model presented in Ref. and, at the same time, provide the simplest setting for studying 
the collision of a kink with a boundary. 
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Time and spatial coordinate are denoted by t and x, respectively, where the latter is taken on a 
circle, that is, in the range x £ [-L, L] (with the endpoints of the interval identified) subject to the Z2 
orbifolding generated by x —* —x. The resulting one-dimensional orbifold S 1 /Z2 has two fixed points 
at x = x\ = and x = X2 = ±L and points x and —x are identified. Instead of working with the full 
orbicircle, we, therefore, can (and frequently will) use the "downstairs" picture where x is restricted 
to the interval x G [0,L]. The fixed points at X{, where i = 1,2, can now be thought of as spatial 
boundaries. 

The field content of our model consists of a single real scalar field 4> with a potential V = V(cp). 
Gravity and gravity-like fields will not be considered here since they are not essential to analyse the 
kink-boundary collision. We require the potential to have at least two different minima (with vanishing 
potential values) in order for kink solutions to exist. We also define the "superpotential" W by 

W' = y / 2V{(/)) . (2.1) 

where here and in the following the prime denotes the derivative with respect to 4>. We will be interested 
in actions of the type 

S = J dt |* dx I \ (^j 2 - \ (^j 2 - V{cp) - 25{x - 0) W{4>) + 26(x - L) W{<j>) J (2.2) 
The equations of motion and boundary conditions in the boundary picture are then given by 



d 2 6 d 2 - 



W'(cf>)\ . (2.3) 



dt 2 dx 2 ^ dx 
For the numerical simulations and for concreteness we will often use the quartic potential 

V(4>)= 1 -g 2 (^-v 2 ) 2 (2.4) 

with associated superpotential 

W = -lg(l<j> 3 -<t>v 2 ) . (2.5) 



2" \3 

The form of the boundary terms in the action (|2.2|) has been chosen to facilitate the existence of BPS 
solutions. This can be understood by focusing on static configurations <j) = 4>(x). For such configurations, 
the action can be written as 

s ~£H%- w ') 2 ■ (2 ' 6) 

Note that the boundaries terms which arise from partial integration when converting the action (|2,2j) 
into the above form are precisely cancelled by the original boundary terms in (|2.2|h Hence, we find 
static solutions cj) = <p{x) of our action satisfy the first order equation 



Conversely, solutions to this first order equation clearly satisfy the boundary conditions as well as the 
second-order equation of motion in (|2.3|) . the latter by virtue of the relation between V and W . 
From eq. (|2.7|1 . the solution </>k of a kink can in general be written as 

<j> K { X )=f{x-Z), f- 1 (</,): 



h W'{4>) 

Here we take 4>o to be the <f> value at the core of the kink, corresponding to the maximum of the potential. 
Then, the constant Z can be interpreted as the position of the kink's core as long as Z £ [0,L]. For 
Z outside this range, the core of the kink is no longer within the physical part of space and Z merely 
indicates the "virtual" position of the core. 

For the quartic potential Q2.4JI the kink solution is given by 

<fna{x) = v tanh (^f) (2-9) 



where the width Lk of the kink is given by 



L K = — . (2.10) 

gv 



As long as kinks are sufficiently away from the boundaries, the field <f> will be close to minima of 
its potential V on the boundaries. Then, the boundary superpotential takes values from a discrete 
set consisting of the values W assumes at the minima of the potential V. These discrete values can 
be interpreted as boundary charges. They are analogous to the topological boundary charges which 
appear in the related M-theory model [21 [I]- I n order to consistently maintain this interpretation, the 
superpotential must be monotonically increasing, and so Eq. (|2.5|) is only valid for —v < (ft < v. 

We note that our boundary conditions are related to the family considered in jUJ. These papers 
studied integrable (1 + l)-dimensional field theories on the half line, including the sine-Gordon theory, 
with a two-parameter set of boundary conditions eft' = MW'((f> — 4>o) at the origin. The fact that these 
boundary conditions do not spoil the integr ability generated considerable interest. 



3 Moduli space evolution 

In the previous section, we have presented a one-parameter family of static BPS kink solutions for our 
theory (|2.2|) which is labelled by the parameter Z. Solutions exist for all values Z £ [—00,00] but only 
for Z 6 [0, L\ does Z correspond to the position of the kink's core. Otherwise, the kink's core is located 
"outside" and only its tail remains within the physical part of space. Values Z = or Z = L correspond 
to a kink with its core being exactly located at one of the boundaries. 

We would now like to consider physical motion in this moduli space of kinks by promoting Z to a 
slowly- varying function of time, Z = Z(t). A time evolution Z — > 0+ (Z —* L_) then describes the 
collision of the kink with the boundary at x = x\ = (x = X2 = L) which is the process we would like 
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to study in this paper. For definiteness, we will focus on collisions with the boundary at x = in the 
following but the conclusions, of course, apply to collisions with the other boundary at x = L as well. 

The one-dimensional effective theory S e s for the modulus Z = Z(t) can be obtained by inserting 
the kink-solution (|2.8[) into the action (|2.2|) and integrating over the spatial coordinate. This results in 

5 eff = J dtF(Z)Z 2 , (3.1) 

where the function F = F(Z) is given by 

F(Z) = W(f(L-Z))-W(f(-Z)) . (3.2) 

The last equation means that F is the difference of the boundary superpotentials evaluated on the kink 
solution. From this interpretation the qualitative behaviour of F can be easily inferred. For Z 6 [0, L] 
and well away from the boundaries (that is, away by more than the width of the kink) the field 4>k is 
very close to neighbouring minima of the potential on the two boundaries. More precisely, (/)k( x = 0) 
is very close to the minimum with the smaller superpotential value and 4>k( x = L) is very close to 
the neighbouring minimum with the larger superpotential value. This means, for such values of Z the 
function F(Z) is finite positive and approximately constant. On the other hand, for Z < or Z > L 
and away from the boundary, <j> is very close to the same minimum for both boundaries which means 
that F(Z) is approximately zero (although still positive). 

The function F can be computed for any given potential V and associated superpotential W to 
confirm this expectation. For the quartic potential 1)2.4(1 . the associated superpotential (|2,5|) and kink 
solution (|2.9|) one finds 

F(Z) = gv 3 

This function has been plotted in Fig. ^ an d it shows indeed the general properties discussed above. 
The constant value of F for Z inside the interval is given by 

F(Z) ~ 4gv 3 /3 for Z » L K and L - Z > L K . (3.4) 

For later purposes, it is also useful to note the asymptotic behaviour of F for Z <C — which will be 
relevant after a collision of the kink with the boundary at x = 0. It is given by 

F{Z) ~ Agv 3 (1 - e~ 29vL ) e 2gvZ for Z <C -L K . (3.5) 

If the width of the kink Lk is much smaller than the size L of the interval, Lk <C L, this reduces to 

F(Z) ~ Agv 3 e 2avZ for Z < -L K and L K <C L . (3.6) 



tanh(— v (x 



Z))- -tanh 3 (^-v(x- Z)) 



(3.3) 



The equations of motion derived from the effective action 1)3.1(1 are simply 



' FZ 2 ) = , (3.7) 



dt 




2 0.2 0.4 0.6 0.8 1 1.2 



Figure 1: The function F(Z) for an interval with length L = 1 and a kink width of Lk = 32 ^ . 



where the dot denotes the derivative with respect to time. This leads to the first integral 



FZ 2 = u 2 = const , 



(3.8) 



where u is a constant. This relation, of course, represents energy conservation. Together with the 
properties of F stated above it implies that Z evolves with constant velocity in the interval Z £ [0, L] 
and, once it has collided with the boundary at x = 0, accelerates and evolves towards Z —* — oo while 
Z diverges. Moreover, from the explicit form of F one can show that Z reaches — oo in a finite time. Of 
course, the validity of the effective theory Q3.1JI when Z and Z diverge after a collision is at best doubtful. 
Although the effective theory provides a good description of slow kink motion for finite values of Z, it 
seems difficult, therefore, to draw any reliable conclusion about the outcome of the collision based on 
this effective theory. Three possibilities are conceivable. The kink could be absorbed by the boundary 
with its kinetic energy being completely converted into radiation (that is, fluctuations of <p around one 
of its vacuum states). In this case the charge on the absorbing boundary would have changed by one 
unit between the initial and final state. It could be inelastically reflected with part of its kinetic energy 
being transferred to radiation or it could be elastically reflected. Then the boundary charge remains 
unchanged. We will now analyse the full two-dimensional theory (|2.2|) to determine which of these three 
possibilities is actually realised. 
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4 Perturbation spectrum 



It is instructive to examine how the eigenfunctions and eigenvalues of small perturbations around the 
kink change with the modulus Z . To find the eigenvalues we must expand the field around the kink 
configuration <p{x) = 4>k( x — Z) + (f(x) and solve the eigenvalue problem 



Here, we have neglected the contribution from the boundary at x = L which is justified as long as 
the kink is sufficiently far away from this boundary. This can always be achieved by making L large. 
We would then expect that at large positive Z, the spectrum would be asymptotically close to that 
of an isolated kink |12j . For the quartic potential (|2.4|) one expects one zero mode with no nodes, a 
bound state with one node and eigenvalue to 2 = 3fi 2 /2, and a continuum starting at u)\ = 2/i 2 , where 
[i = gv/\/~2. At large negative Z, the spectrum should be asymptotically close to that of the vacuum, 
which also possesses a nodeless zero mode and a continuum starting at uj 2 = 2/u 2 , but no bound state. 

Although the free kink eigenfunctions and eigenvalues are known exactly |12j , the non-trivial bound- 
ary condition, or equivalently the extra <5-function contribution to the potential, change the problem. 
As we have not been able to find a solution to the eigenvalue problem in closed form, we resorted to a 
simple numerical method. 

Eq. (|4.1I) was discretised on an 2N + 1-point interval — L < x < L, using the lowest order discrete 
Laplacian 



where Ax = L/N. At the left (right) boundary the double forward (backward) derivative was taken, 
which improves the accuracy at the expense of generating two spurious eigenfunctions with large negative 
eigenvalues. The <5-function was defined as 5{x) — > 5io/Ax. One should note that in defining the problem 
for both positive and negative values of x, we will find eigenfunctions which are both symmetric and 
antisymmetric under reflections around x = 0. We are clearly interested only in the symmetric ones. 

The eigenfunctions and eigenvalues are easily found using MATLAB. We took N = 100 and L = 10, 
in units where fj, = 1. The lowest three (0, 1 and 2 nodes) are shown for —2 < Z < 2 in Fig. @. The 
closeness of the lowest eigenvalue to is a measure of the accuracy of the procedure. The shapes of the 
eigenfunctions for Z = —2, Z = and Z = 2 are shown in Fig. (J3J). 

As we will see later, the existence of a zero mode for all values of the kink position Z is crucial for 
an understanding of the kink-boundary collision. We will also find that it has important implications 
for the creation of kinks due to spontaneous emission from the boundary. 






Ax 2 



(4.2) 
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Kink position Z 



Figure 2: The lowest three eigenvalues uj 2 as a function of kink position Z. As Z decreases, the bound 
state at uj 2 = 1.5 moves towards the continuum at u> 2 = 2. In infinite volume it will remain bound, but 
exponentially weakly. 
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Figure 3: Graphs of the eigenfunctions corresponding to the lowest three eigenvalues for the kink at (a) 
Z = —2, (b) Z = 0, and (c) Z = 2. At the top is a graph of the field cf> itself. 



5 Numerical Simulation 



As we have seen in the previous sections, a straightforward moduli calculation suggests that, as the 
kink collides with the boundary, it should cross it and accelerate towards infinity. As a result, the 
field configuration inside the bulk will approach one of the true vaccua of the theory. We could then 
expect the kinetic energy of the kink to be converted into radiative excitations of the vacuum that 
would propagate away from the boundary. This reasoning suggests, the kink should be absorbed by the 
boundary as a result of the collision. In this section we will test this expectation by performing numerical 
simulations of the kink/boundary collision. As we will see, a more complex picture will emerge, showing 
how the moduli description fails to capture the essential features of the collision process. 

We model the bulk equation of motion Eq. Q2.3JI with the quartic potential (|2.4|) on a ID lattice 
using a second order discretization of the Laplacian operator, and propagate it in time according to 
a standard leapfrog algorithm. The implementation of the boundary condition is less straightforward, 
since it relates a function of the field to its spatial derivative. Taking as an example the boundary 
condition at x = 0, we discretise it as 

= W'(0 m ), (5.1) 



8x 

where <j) m is the value of the field evaluated at some point in the interval (0, 5x). We allowed for 
different choices of the discretization scheme by writing <f) m = (1 — G) <j)Q + G<j>\, with G a parameter 
lying between and 1. Since the discretised spatial derivative is effectively evaluated in the mid-point 
between lattice sites i = and i = 1, the most natural option is to set G = .5 We confirmed that this 
choice is the most accurate one, by evolving in time a stationary kink configuration, with its core placed 
near the boundary. For G > .5 the kink is slightly repelled by the boundary, whereas for G < .5 the 
boundary has an attractive effect. The disadvantage of setting G = .5 is that it turns Eq. (|5.1|) into 
a non-linear implicit equation for 0o, which has to be solved numerically for each time-step. This is 
clearly the case for any choice of G other than G = 1. On the other hand, for the lattice action obtained 
by a straightforward discretization of Eq. (|2.2j) . the only energy conserving algorithm corresponds to 
G = 0. For G > we cannot rely on energy conservation as test of the adequacy of the choice of 
time-step. Of course all schemes conserve energy in the 5x — > limit. Since we are dealing with a one- 
dimensional system, we had no serious computational restrictions and we chose to work with G = .5. 
We have checked that, in the limit of small lattice spacing the three values G = 0, .5 and 1 lead to the 
same results. The boundary condition equation Eq. (|5,1|) was solved using an implementation of the 
Newton-Raphson method (see \T3\). 

The initial conditions for the numerical evolution should model, as closely as possible, a kink in 
the bulk with a given velocity Zi moving towards one of the boundaries. We set up the initial field 
configuration as <Pk{x — Zi) using the kink profile in Eq. ()2.9j) . with Zi the initial core position far from 
both boundaries. The initial field momentum is defined as II = — Zi d x (f>K, where Zi is the chosen initial 



kink velocity. The simulations were run for a theory with g = 2y/2 and v = 1, which implies a kink 
width of Lk = We chose a box of physical length L = 16, comfortably larger than the kink width. 
For high and mid-range kink velocities, the lattice parameters were set to 5x = .0125 and 5t = .005. 
Low collision velocities lead to a large total simulation time and a smaller time-step is needed to avoid 
error accumulation. For the lowest velocities we used 5x = .00125 and 5t = .001. 

We explored a wide range of collision velocities, letting Zi vary between .01 and .75 in steps of .01 
In Fig. 0]we can see a series of snapshots of eft for several evolution times in a typical run (in this case 
Z{ = .3). During the first part of the evolution the kink moves towards the boundary with constant 
velocity. As its core starts moving into the boundary there is, as expected from the moduli picture, 
an increase in Z. At some point the kink leaves the bulk completely and the field is left in the true 
vacuum, oscillating rapidly in the vicinity of the boundary. The kink then "returns" into the simulation 
box, moving away from the boundary with decreasing velocity at first. When its core is well inside the 
bulk the kink velocity oscillates around an average value Zt. This final value is always smaller than the 
initial velocity Zi, as some energy is lost into radiation during the collision. In the late time profiles in 
Fig. 0]we can see these radiative perturbations moving away from the boundary. 

In Fig. 0we have a plot of the kink velocity versus time. Z(t) was evaluated by numerical differ- 
entiation of the position of the core of the kink Z(t). This in turn was defined as the position of the 
zero of the field, obtained by interpolating <f> as its sign changes. The same features of the evolution can 
be seen, with Z(t) increasing as the kink approaches the boundary, and being reflected later on with a 
lower final velocity. Note that in both the beginning and end periods of the evolution, the velocity is 
not strictly constant, oscillating around an average value. This indicates that in both cases, apart from 
the translational mode, other modes of the kink are excited. As far as the initial condition is concerned 
this was expected since we are not using a Lorentz boosted kink, which would clearly not obey the 
boundary conditions. This effect is particularly visible for high kink velocities, both before and after 
the collision. 

Surprisingly, the overall qualitative pattern of the evolution did not change within the whole range 
of velocities observed. Even for very low velocities, where we expected the reasoning based on the 
moduli calculation to be applicable, the incoming kink was always reflected by the boundary after the 
collision. In the next few sections we will try to explain this behaviour and gain a deeper quantitative 
understanding of the collision process. 

6 Collision as a vacuum perturbation 

As we have seen in Section |31 the motion of the kink should be well described by the moduli equations 
of motion, up to the point when the velocity of its core becomes too large. For low values of the initial 
velocity Zi, we expect this to happen when the kink core is outside the bulk and far from away from 
the boundary. That this is the case can be easily shown using the moduli energy conservation equation 
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Figure 4: Kink collision for initial velocity Zi = 0.3. In the first two snapshots, t = .2 and t = 10, 
the kink moves towards the boundary with approximately constant velocity. For t = 18 most of the core 
has already crossed the boundary, although the profile inside the bulk remains relatively undistorted. For 
t = 19 the field oscillates around the true vacuum. These perturbations propagate away at the speed of 
light, as seen in the next two plots. At t = 22 the kink core has re-entered the bulk, and in the final plot 
we can see the final kink profile, moving away from the boundary with constant velocity. 
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Figure 5: Kink velocity versus time for initial velocity Zi = 0.3. As predicted by the moduli calculation, 
the kink accelerates as it approaches the boundary. After the collision it returns into the bulk with a 
lower final velocity. During the collision period, when the core is outside the bulk, the kink position 
cannot be determined by the zero of the field, and we have defined the kink velocity as zero. 
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(|3.8|) . Assuming that the moduli approximation breaks down for velocities above a certain Z max . Then 
Eq. (jSIBj) implies that the corresponding core position Z max should satisfy 

F{Z max ) = ~gv 3 (^-] • (6.1) 

Here we assumed the kink to be initially inside the bulk and far from its boundaries so that (|3.4j) applies. 
Clearly, the form of F(Z) implies that for Zj small enough, Z max can be made arbitrarily large. 

If, when the maximal velocity is reached, we have Z max <C — £k> the field configuration in the bulk 
will be everywhere very close to the vacuum v. This suggests that for low Zj, we could use the linearised 
equations of motion around v to describe the evolution, after the moduli approximation breaks down. 
In this section we will use this approach to show that the kink bounces back inelastically after the 
collision, confirming the previous simulation results. 

6.1 Expansion basis for linearised equations of motion 

As discussed above we perturb around the vacuum v and write 4> = ip + v. For small 92 the field equation 
of motion reduces to 

(p = tp" - m 2 Lp (6.2) 

where m = gv in the case of the quartic potential. The non-trivial contribution to the evolution will 
come from the boundary condition, which is given by 

<p'{xi,t) = (g/2) I 4> 2 {xi,t) - v 2 \ = m\(p(xi,t)\ (6.3) 

In this section we will focus on a collision taking place at the left boundary, x = 0, and effectively 
assume the bulk to be half-infinite. This should be a reasonable approximation as long as the size of 
the bulk is larger than the collision time, so that effects from the radiation produced during collision 
reaching the opposite boundary can be neglected. Note that, as a consequence of the norm in ()6.3j) . the 
boundary condition for </? is non-linear. Nevertheless, for periods of the evolution during which 93 does 
not change sign at the origin, (J6.3|) reduces to a linear equation. In particular we have 

c^'(0,t) = m<p(0,t), if ip(0,t)>0, (6.4) 
tp'(Q,t) = -mip(Q,t), if tp(Q,t) <0. (6.5) 

This implies that we can still expand arbitrary solutions of (|6,2j) and ()6.3|) in terms of given mode 
functions, as long as we use a different basis for different periods of the evolution. We start with the 
case 99(0, t) > 0. A complete basis of solutions of ()6.2j) satisfying (J6.4|) is given by 

ip k (x, t) = N k e ±iult [cos(fcr) + — sin(A;a;)j , w 2 = k 2 + m 2 , k> . (6.6) 
L k J 

To obtain the normalisation coefficient we define the mode functions in a finite box of size L, leading 
to iV~| = 2k 2 /(Luj 2 ). In the end of the calculation we will take L — > 00. As for k, it will have a discrete 

12 



spectrum depending on the boundary conditions on the right end of the box. For simplicity we chose 
periodic boundary conditions, implying cos(/cL) + (m/k) sin(kL) = 1. As the box size is taken to infinity 
this should not be relevant for the final result. 

In the case of negative field values <£>(0, t) the basis is given by 



m 



ip k (x, t) = N k e ±tujt cos(kx) - — smikx) , w 2 = k 2 + m 2 , k>0 (6.7) 
L k J 

with cos(kL) — (m/k) sin(kL) = 1. Here N? is a normalisation which will be unimportant for the 
subsequent calculation. In addition, for the case of negative field values cp(0, t) and only for this case, 
there is also a zero-mass solution of the general form 

ip (x,t) = {a + bt)e- mx (6.8) 

This zero mode corresponds to the tail of the kink when its core is far outside the bulk. It can be seen 
as the linearised version of the translational zero mode of the kink. As we will see, this is the essential 
feature that distinguishes the evolution in the two regimes, ip > and ip < 0. 

6.2 Linear solution 

We now apply the above results to the kink/boundary collision situation. For simplicity, we define 
t = as the time when the moduli approximation breaks down. For this time, we must find the 
corresponding field configuration inside the bulk which can then be taken as initial condition for the 
equations of motion linearised around the vacuum. When the kink core is very far from the boundary, 
the bulk field and its time derivative can be approximated by 

4{x,t) = v - 2v e mZ ® e~ mx (6.9) 
d t <t>{x,t) = -2mvZ(t)e mZ(t) e~ mx (6.10) 

We will assume that we can neglect the small perturbation around the vacuum and that effectively 
<j) = v in the bulk at t = 0. The derivative term on the other hand, remains finite even as Z(t) diverges. 
Using energy conservation (|3.8|) and (|3.6f) we obtain 

Amv 2 e 2mZ Z 2 = F{Zi) Z* (6.11) 

which implies, using the asymptotic expansion for F(Z), that 

e mZ i=y. (6.12) 
v3 

This allows us to express the finite amplitude of the velocity profile at collision time in terms of the 
initial velocity of the kink. The perturbation field at t = is then defined as 

<p(x,0) = (6.13) 

<p(x,Q) = — ^mvZ ie - mx (6.14) 
v3 
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These are now to be taken as initial conditions of the linearised equations of motion. Since ip > we 
will be entering a phase of the evolution for which <p(0, t) > and we should expand (|6.14j) in terms of 
the corresponding basis (|6.6|) , The expansion coefficients are easily obtained (the result depends only 
on integrals of products of exponentials and trigonometric functions) leading to 

ip(x,t) = ——m (l — e~ mL ) — =■ sin(iL>t) cos(A;a;) + — sin(A;x) (6.15) 
L ^— ' uj° V k 

k 

where A = (2/y/S) mv Zi is the amplitude of the exponential tail in the initial velocity profile for the 
field. 

Equation (|6.15j) is only valid as long as ip(0, t) > 0. When the field becomes negative at the origin, 
the boundary condition changes sign and the solution should then be expanded in terms of the second 
basis (|6.7j) . Let us define a "reflection time" tR, for which the solution becomes negative at the boundary, 
that is (p(0,tji) = 0. We can use (|6.15|) to calculate explicitly the mode amplitudes of ip(x,tji) = in 
terms of the <p> < basis. The new ingredient now is the presence of the zero mode. All other modes 
are oscillatory, corresponding to the radiation that will propagate away from the boundary after the 
collision. The zero mode tpo on the other hand, will grow linearly in time, becoming the tail of the 
incoming kink. The situation will now mirror what happened before the collision, with the amplitude 
of the massless mode being related to the final velocity of the reflected kink Zf. We define the zero 
mode amplitude as 

a = J dx dt<p(x, t R ) N e~ mx (6.16) 

where Nq = 2m/(l — e~ 2mL ) is the corresponding normalisation factor. The time derivative of the field 
at tR will be given by: 

(p(x, tfi) = ao Nq e ~ mx + oscillatory terms (6-17) 

Since we are mainly interested in determining the final velocity of the reflected kink, we will not calculate 
explicitly the positive k component of the field. We define the amplitude of the exponential term in 
ip(x, tfi) as B = ao No, and in analogy with the calculation for the kink before the collision, (|6.14l) leads 
to 

Z f = ^-B. (6.18) 

A reflection coefficient can be defined as R = Zf/Zi = B/A, quantifying the elasticity of the colli- 
sion. The amplitude B can be easily obtained by performing the integral in Eq. H6.16|) . leading to the 
expression 

H=7 = J»?J(1- e- mL ) 2 4m 2 £ ^ 008(0**) (6.19) 



A L ' ^ u G 

k 

for the reflection coefficient where u is given in Eq. (|6.6j) . Note that this solves the problem. We have 
obtained a formula for the reflection coefficient in terms of the reflection time tR, which, in turn, is 
implicitly determined by Eq. Q6.15|) as the time satisfying ip(0, tR) = 0. 
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6.3 Comparison with numerical data 

To evaluate exact numerical values it is convenient to approximate the series obtained in the previous 
section by integrals. There is a slight subtlety involved, related to the fact that we must add over all 
values of k that solve the periodic boundary condition cos(fcL) + (m/k) sin(/cL) = 1. This equation 
is solved by k = 2nir/L, and by an extra non-trivial solution in every interval (2nir/L,2(n + 1)tt/L). 
This second solution approaches k = 2nir/L as n — > oo. In replacing sums over k by integrals we 
must take into consideration the fact that there are two k's in each interval of length 2n/L and use 
(vr/L) J2k J dk (rather than (2ir/L) — > J dk). Using the continuum version of (|6.19|) we obtain 
for the reflection coefficient 

16m 3 k 2 

R = / dk — r cos(ujt R ) (6.20) 

vr Jo w 

The reflection time tR is given implicitly by 

r°° k 2 

J dk—^ sxu{ut R ) = (6.21) 
Both integrals can easily be re-scaled and evaluated numerically, leading to the final results: 

= __ R = o 63 , 622 \ 

m 

These compare very well with the simulation results in the limit of low collision velocities, as illustrated 
in Figs. and [7| As the initial kink velocity increases deviations from the analytical result can be 
observed. This is to be expected, since for higher velocities the moduli approximation should break 
down earlier, making equation (|6.14|) a worse guess for the form of the initial conditions for the linear 
regime. Also, we should expect corrections to the linear approximation as the field starts probing the 
non-linearities of the potential during the collision, for high Zj. Still, in the particular system simulated, 
the numerical results remain within 10% of the analytical prediction for velocities up to Zi = 0.3. 

It is interesting to note that the equation 1)6.22(1 depends only on the mass of the perturbed theory 
around the vacuum. In fact, it is easy to see that the derivation leading to this result could be reproduced 
for a general potential with non-vanishing second order derivative at the true vacuum. In particular, 
for this class of potentials, the reflection coefficient should be model independent. We confirmed this 
remarkable result by simulating a few collision velocities for a sine-Gordon model. As expected, we 
observed R = 0.63. 

7 Moduli space interpretation of the kink-boundary collision 

Can we find an interpretation for our results in terms of a suitable modification of the one-dimensional 
moduli-space approximation described in Sectionl^P As before, we focus on a collision with the boundary 
at x = 0. In view of the effective action (|3.1|) . it seems natural to introduce a new field C with a 
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Figure 6: Reflection time versus collision velocity for a model with g = 2\/2, v = 1, m = 2\[2. The 
theoretical "prediction \6.2 6 2\) for this choice of parameters is tn = 0.73. The reflection time tn was 
defined as the length of time for which </>(0, t) > v for each simulation. 



canonically normalised kinetic term by setting 



dz 



±VF 



dZJF(Z) 



(7.1) 



Note that the integral in this definition is indeed finite which corresponds to the fact that Z reaches 
— oo in a finite time. For a positive (negative) sign in the above definition the full range of Z values is 
mapped into positive (negative) ( values while the point Z = — oo is mapped to ( = for both signs. 
The effective action 1)3. 1JI now takes the canonical form 



S, 



eff 



/ 



'-2 



and leads to the first integral 



( = V = const . 



(7.2) 



(7.3) 



Let us discuss the collision process in terms of the new field £. We start off at some finite, positive value 
(, corresponding to the kink being located between the boundaries, and move with a constant, negative 
velocity towards £ = 0. When £ reaches some smaller, positive value, corresponding to the position of 
the boundary at x = 0, the collision takes place. From there £ continues to evolve to zero and, finally, 
into negative values. Nothing special happens at £ = which corresponds to the formerly problematic 
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Figure 7: Reflection coefficient versus collision velocity for g = 2y/2, v = 1, m = 2y/2. The (model 
independent) analytical prediction is R = 0.63. Both the initial and final velocities were determined by 
fitting the position of the core of the kink ( as defined in Section [TJJj to a straight line in the relevant 
regimes. 
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point Z = — oo. We would now like to re-interpret this evolution of £ in terms of our original field Z. 
Note that, from Eq. Q7.1|) we have 



Z = 



+7F fOT C>0 



(7.4) 



This means crossing of £ = after a collision implies a sign change in the velocity of Z and, hence, 
a reflection of the kink. This is precisely what we have found in the full two-dimensional theory. An 
alternative way of stating this is to note that £ is really defined on R/Z2. Here the Z2 identifies £ and 
— C which correspond to the same field value Z. In fact, upon including the other boundary the £ moduli 
space is, perhaps not surprisingly, given by S 1 /Z2- Note, however, the fix points of this moduli space 
orbifold do not quite correspond to the fix points of the space-time orbifold (which are at x = 0, L) but 
rather to the points "at infinity", that is, Z = ±00. 

There are two obvious conclusions about the collision process which we can deduce from our simple 
picture. First, the absolute values of the kink velocity, |V|, before and after the collision are the same, 
that is, the reflection is elastic and the reflection coefficient is R = 1. This is in obvious disagreement 
with our two-dimensional results which showed the reflection is inelastic even at low velocities with a 
reflection coefficient R ~ 0.63. This disagreement arises because the modified moduli space picture still 
does not correctly describe the system for a short time period, corresponding to the reflection time tn, 
when the modulus £ is not defined. During this short period, the additional radiation is created, a 
process which clearly cannot be described by our effective moduli space theory with a single degree of 
freedom. 

Secondly, the total time t ie Q between the two crossings of the core of the kink with the boundary is 
given by 



Note, this time is different from the reflection time tn defined earlier and, for low velocities, we have 
^refl 3* tR. While i re fl is the total time the kink vanishes "behind" the boundary, tn is the time the 
moduli space approximation breaks down when the kink is far outside at Z — > — 00 or £ ~ 0. Comparison 
with the numerical simulations shows that the above formula approximates t re a well for low velocities 
with the typical discrepancy being of order tR. This is understandable since the main contribution to 
the integral in Eq. (|7.5jl comes from the region Z ~ around the boundary where the moduli space 
approximation is still valid. 

Another conclusion from Eq. (|7.5jl is that the kink can disappear behind the boundary for an 
arbitrarily long time for sufficiently low initial velocity. However, as we have seen, it eventually always 
reemerges and is reflected back into the bulk. 








(7.5) 
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8 Conclusions 



Combining the numerical and analytical results discussed above, a consistent picture of the kink/boundary 
collision process begins to emerge. We found that for most of the evolution the system is well described 
by an effective theory for the kink position modulus. As the kink leaves the bulk and the modulus 
diverges, there is a short period of time (the reflection time £r, in the notation of Section EJ) for which 
this approximation breaks down. During this stage, the field in the bulk takes values outside the interval 
(—v, v) between the two minima of the potential and, as a consequence, the moduli approximation is not 
defined. Solving the full field theory for this period of the evolution, we showed that the kink is always 
inelastically reflected back into to the bulk. After the collision the kink moves in the bulk following once 
again the moduli equations of motion. We presented an analytical approach which describes the kink 
reflection as a fluctuation of the scalar field around its vacuum state. This picture led to analytic expres- 
sions for the reflection coefficient R (the ratio of initial and final kink velocity) and the reflection time 
tfi. These analytical results turned out to be in good agreement with the numerical simulations for low 
kink velocities, thereby supporting our analytical approach. Remarkably, the reflection coefficient and, 
hence, the fraction of energy lost into radiation during a collision, turned out to be model-independently 
given by R ~ 0.63. 

Although we have shown that the reflection of the kink can be modelled by a formal continuation 
of the moduli space, some aspects of the collision, such as its inelasticity, and an extra time delay 
tn, are not correctly reproduced within this moduli space approach. So far we have not been able to 
find an alternative description of the collision that does not require taking into account the full set of 
degrees of freedom of the two-dimensional field theory. It is possible that a more adequate moduli space 
description can be achieved by a different continuation of the moduli space but clearly more work is 
needed to settle this question. 

These results may be disappointing at first as they seem to imply that the number of kinks in the 
bulk and the charges on the boundaries are always conserved during collisions. Initially we had expected 
to find examples of collisions for which the kink would be absorbed by the boundary, leading to a change 
of the boundary charge pQ. In the related M- theory model such a change in the boundary charge 
indicates a change in the topology of the M-theory background and is, therefore, of particular interest. 

However, on closer examination, the mechanism behind the reflection of the kink reveals a much 
richer potential for processes changing the boundary charges. As shown in Section H3 after the collision 
the kink grows back into the bulk as the exponential zero mode around the vacuum becomes excited. 
If we now consider a general arbitrary perturbation of the vacuum, it is very likely that it will have 
a non-vanishing component in the direction of the zero mode. This will lead to a kink forming and 
moving into the core, changing by —1 the charge of the boundary. The somewhat unexpected conclusion 
seems to be that it is reasonably easy to extract a kink from the boundaries. In fact it is easy to show 
that this process is not inhibited by a lower energy threshold. When a kink forms, the change in the 
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bulk energy corresponding to its core is compensated for by an equal change in the boundary energy 
terms. The energy of the initial perturbation is then converted exclusively into kinetic energy of the 
incoming kink. This implies that kinks with arbitrarily low velocities can always be obtained from field 
configurations arbitrarily near the vacuum. In this sense emission of kinks from the boundaries which 
change the boundary charges becomes unavoidable, in practise. 

Let us illustrate these ideas by looking at the very long time limit of the kink/boundary collision 
process. Here we will follow and consider a theory with a periodic potential V{4>). In the situations 
considered so far the field values were always near the two vaccua —v and v, and the large 4> structure 
of the potential was never relevant. This is not the case, as we will see, for considerably long times after 
the collision. As we have discussed, when the kink collides with the left boundary some radiation is 
produced which propagates away in the direction of the opposite right boundary (and vice versa). As 
the vacuum in the right boundary becomes perturbed, we can expect as before, that the corresponding 
zero mode becomes excited. From our arguments above this zero mode should with time grow into 
an additional, left-moving kink. Something similar could be expected to happen in the left boundary, 
after the reflected kink moves away into the bulk. Since we expect some small perturbation around 
the vacuum there, another (probably very slow) right-moving kink should start forming. For a periodic 
potential and very long times this process would repeat itself several times, with increasing numbers of 
kinks being extracted from both boundaries. That this seems to be indeed the case has been confirmed 
in Fig. |HJ For a sine-Gordon model with potential 



and parameters m = 1, A = 1 we observed the formation of five kinks from an initial single-kink collision. 
As expected, the kinetic energy of the initial kink is distributed amongst the final kinks that move with 
much lower velocities. 

As an extra test we evolved a number of near vacuum configurations with zero bulk charge and 
observed the formation of a number of kinks at both boundaries, as expected. The natural question 
at this point is which mechanism would cause the kink formation process to stop eventually. One 
possibility is that, as their velocities become increasingly small, a kink gas would form, its density being 
limited by the repulsive kink-kink interaction. We can conjecture that for a very large family of initial 
conditions a gas of stacked branes should form, with density closely related to the energy of the initial 
field configuration. If this happens to be the case, the final charge configuration, in particular the final 
boundary charges of the system, would be selected by its total energy. Clearly a deeper understanding, 
both in numerical and analytical terms, of the system is needed in order to clarify these intriguing 
possibilities. 
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Figure 8: Initial single colliding kink gives rise to configurations with increasingly higher kink number. 
A sine-Gordon potential (|8.1|) with m = 1 and A = 1 was used. The y-axis is normalised such that the 
vacuum points correspond to integers. The initial kinetic energy of the single kink is distributed amongst 
the final kinks that have lower velocity. 
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